Stochastic Gravitational Wave Background originating from Halo Mergers 
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The stochastic gravitational wave background (GWB) from halo mergers is investigated by a 
quasi-analytic method. The method we employ consists of two steps. The first step is to construct 
a merger tree by using the Extended Press-Schechter formalism or the Sheth & Tormen formalism, 
with Monte-Carlo realizations. This merger tree provides evolution of halo masses. From A'^-body 
simulation of two-halo mergers, we can estimate the amount of gravitational wave emission induced 
by the individual merger process. Therefore the second step is to combine this gravitaional wave 



emission to the merger tree and obtain the amplitude of GWB. We find flaw 



10" 



for / 



10" 



10" 



Hz, where Qaw is the energy density of the GWB. It turns out that most of the 



contribution on the GWB comes from halos with masses below 10^^ Mq and mergers at low redshift, 
i.e., 0<z < 0.8. 
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I. INTRODUCTION 



Stochastic gravitational wave background (GWB) is a 
valuable source of information for both cosmology and 
astrophysics. 

In the very early universe, inflation [iHJ is expected 
to be the most feasible scenario. It is believed that both 
scalar and tensor mode perturbations from quantum fluc- 
tuations are produced in the epoch of inflation. The for- 
mer evolve into large scale structure of the universe. And 
the latter directly travel toward us and are observed as 
the GWB. Therefore this GWB can be a direct probe 
of inflation. For that, one needs to solve the transfer of 
the gravitational wave in the history of the expanding 
universe from the epoch of inflation to present. For ex- 
ample, Kuroyanagi et al. Q studied evolution of the grav- 
itational waves from inflation in detail and showed that 
the amplitude of the GWB energy spectrum is ^ 10^^^ 
for / > 10^^ 'Hz while the value is very dependent on 
models of inflation. 

The GWB is also produced by second-order scalar per- 
turbations. In the linear perturbations, it is known that 
structure formation does not generate any gravitaional 
waves. It turns out that, however, gravitational waves 
with cosmological scales would be emitted during the 
process of structure formation if one includes the second- 
order terms in the perturbation equations as is studied in 
. The result obtained by them is that the amplitude 
of the GWB is 10~^° — 10~^^ for a wide frequency range. 

The best method for detecting such GWB with cos- 
mological scales is to measure the B-mode of the Cos- 
mic Microwave Background (CMB). More specifically, 
the CMB photons are very sensitive to gravitational 
waves with / ^ 10~^''llz Q. In actuality, several 
missions suited for this aim have been planned; for 



example, ACT Pol, SPTpol, PQL ARBeaR, LiteBIRD 
(sec, http://cmbpol.kek.jp/index-e.html), and Cosmic 
Inflation Probe (see, http://www.cfa.harvard.edu/cip/). 
Their primary purpose is to detect primordial gravita- 
tional waves generated during inflation. 

In this paper, we investigate the GWB originating from 
dark halo mergers, which are expected to produce gravi- 
tational waves with cosmological scales. In a hierarchical 
model of structure formation, it is expected that low- 
mass dark halos repeatedly merge with each other, then 
more massive dark halos are formed (see, e.g., [Tol. flit). 
A large amounts of gravitational waves from the process 
would be emitted because the process is a highly non- 
linear event. In (T^ . we studied the gravitational waves 
emitted from a single galaxy merger with N-body simu- 
lations. The peak luminosity and total emitted energy 
were found to reach about 10"^^ erg/sec and 10"'^ erg, 
respectively for a collision of two galaxies with masses 
3.8 X 10^'^ Mq. We also studied the relative contribution 
of the disk, bulge and halo, the effect of initial velosity 
and relative angular momentum. To calculate the GWB, 
we sum the gravitational wave spectrum from a single 
merger over the merger history. 

This paper is organized as follows: in section 2 our 
method for calculating the GWB is described; in section 3 
we show the results; and in section 4 we provide summary 
and discussion. In this paper we consider a spatially- 
flat CDM model with the following cosmological param- 
eters: n,n = 0.275, r^A = 0.725, /i = 0.702,(78 = 0.815 
and Hg ~ 0.963 for the density parameters of matter and 
cosmological constant, the bubble constant, the density 
fluctuation in sphere of 8h~^ Mpc and the spectral index 
of the primordial scalar fluctuations p^ . 
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FIG. 1: The mass function at present by using EPS(solid) and 
ST(dashed). Tlie y-axis is the number density in the unit of 
l/(h-iMpc)^ The X-axis is the halo mass, M. 



II. METHOD 

In this section we describe our quasi-analytic method. 
As the first step, we construct merger history using 
Monte-Carlo simulation based on the method of [ij]. For 
the mass function of dark halos, we consider Extended 
Press-Schechter formalismjEPS) [H, [ill and Sheth & 
Tormcn formahsm (ST) [Hlill, although the former is 
mostly used below. In the EPS formalism which provides 
the conditional mass function, the fraction of the trajec- 
tories in haloes with mass Mi at zi that are in halos with 
mass Mq at zq (A/i < Mq, zq < ^i), is given by, 



fps[Mi,zi\Mo,zo)d(jl 
1 5cs - 5cfl 



27r (erf 



-13/2 



exp 



dal (1) 



where (5c, i and ctj arc the threshold in the Press-Schechter 
formalism at time Zi and the mass variance for mass Mj , 
respectively. 

On the other hand, in case of ST, conditional mass 
function is described by, 
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(2) 



with A = 0.3222, a = 0.707 and p = 0.3. 

By using Eqs. ((!]) and ([2]), the probability that is a par- 
ent halo with mass Mi at zi being in a progenitor halo 
with mass Mq at zq is obtained and the merger history 
can be constructed by performing Monte-Carlo realiza- 
tion according to the conditional mass function. In Fig. 
[T]mass functions by EPS and ST arc shown. One can see 
that the ST mass function has a larger number of massive 
halos. Below we use EPS unless otherwise noted. 



The merger history obtained here provides only evo- 
lution of the halo mass and does not give spatial infor- 
mation such as the relative velocity and the angular mo- 
mentum of halos, which one may consider as important 
ingredients to estimate the gravitaional wave emission. 
Fortunately, as we showed in [l^ , the initial relative ve- 
locity does not affect the gravitational wave emission so 
much. Specifically, the difference in the emitted energies 
is about 20% between models with zero initial relative 
velocity and 220 km/s which is the maximum initial rel- 
ative velocity of two gravitationally-bound halos. On the 
other hand, the relative angular momentum has a rel- 
atively larger effect on the gravitational wave emission. 
A head-on collision emits three times more energy than 
a collision where two halos initially have a circular or- 
bit. In the current paper, we assume head-on merger for 
all the collisions so that our calculation would be over- 
estimation by a factor less than three. 

In constructing the merger history we take the mass 
resolution Mi = lO^'^h^^ Mq, that is, we ignore halos 
with masses less than Mi and only consider merger of 
two halos. The adopted time step of the merger histories 
is a redshift interval ofAz = 0.06(l-|-z), corresponding to 
the dynamical timcscalc of halos that collapse at redshift 
z. 

As the second step, we sum up the gravitational waves 
from dark halo mergers following the merger history ob- 
tained above, using the result of [l^l which provides the 
gravitational-wave spectrum from a single merger as a 
function of the halo masses. This spectral function must 
cover a wide range of masses and mass ratios and in 
we found a scaling relation of the spectrum with 
respect to them. The spectrum of gravitational waves 
from a merger of equal-mass halos can be written as 
EGW,o{f:Mo) where / and Mq are the frequency and 
the fiducial mass. The energy spectrum for a given mass 
M can be decribed by. 



EGw{f,M) 



M 
Mo 




.f,Mo 



where a = 3.4 
unequal masses 



(3) 



The energy spectrum in the case of 
Ml and Af2, can be obtained by re- 
placing M in the Eq. Q with y/M[M^. Here we set 
Mo = 3.8 X 10^^h~^MQ and then the total energy is 
Egw.o = 5 X 10^6 h-^eig. 

Here wc demonstrate how valid the scaling relation Eq. 
is. Fig. [2] shows the energy spectra of the emitted 
gravitational waves in the logarithmic spacing, i.e, fE{f) 
for three equal-mass mergers with M = 3.8 x 10^^^ 3,8 x 
10" and 3.8 X 10^" /i^^Mq. The top panel is the original 
spectra and the bottom panel is the spectra scaled by Eq. 
© to fit with the case with M = 3.8 x lO^^h'^MQ. In 
these cases, the error of the scaling relation is less than 
40%. It should be noticed that the energy spectrum has 
a peak between 10~^^ to 10~^^Hz which corresponds to 
the dynamical time scale of halo merger, ~ IGyr. 

Fig. [3] is a demonstration for mergers of unequal-mass 
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FIG. 2: The energy spectra of the gravitational waves in 
the logarithmic spacing for mergers of equal-mass halos, 
M/h-^MQ = 3.8 X 10^2 (solid hne), 3.8 x lO" (dashed line) 
and 3.8 x 10^" (dotted). Y-axis is the energy spectra of the 
gravitational waves from a merger in the unit of /i~^erg/s. 
Top: the original spectra. Bottom: the spectra scaled by Eq. 
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FIG. 3: The energy spectra of the gravitational waves in 
the logarithmic spacing for mergers of unequal-mass halos. 
The solid line shows the case of a merger of equal-mass halos 
M/h-^MQ = 3.8 X 10^2 shown for comparison. The dashed 
and doted lines show the cases of unequal masses (dashed 
line:M//i"^M0 = 3.8 x lO" and 3.8 x 10^^, dotted line: 
M/R-^Mq = 3.8 X 10^° and 3.8 x 10^^). Y-axis is the en- 
ergy spectra of the gravitational waves from a merger in the 
unit of /i~^erg/s. Top: the original spectra. Bottom: the 
spectra scaled by Eq. ([Sjl. 
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FIG. 4: The GWB spectra for each mass range. The solid, 
dashed and dotted lines are mass range 
10^^ 10" and 10^^ respectively. 



calculated as, 

PGW[ 



I'com 



(4) 



where Vcom is a comoving volume and Ni represents the 
A^i-th merger in the i-th redshift bin. It should be noted 
that, due to the expansion of the universe, the frequency 
and the energy density of the gravitational waves are red- 
shifted as / oc 1/(1 + z) and pcw oc 1/(1 + z)^, respec- 
tively. Therefore, the energy density pcw at = can 
be written as, 



PGwif, z = 0) = 



where / and Vphys show the redshifted frequency and a 
physical volume, respectively. Finally, the density pa- 
rameter of the GWB, i^GW, is defined as. 



^Gwif) = 



1 



PcfiC 



dpGwif,z = 0) 



din/ 



(6) 



where pc,o and c show the critical density of the universe 
at present and the speed of light, respectively. 



halos. Here we fix the higher mass to 3.8 x 10^^h~^ Mq 
and vary the mass ratio as 1 : 1, 1 : 1/10 and 1 : 1/100. 
The top panel in Fig. [3]shows the original spectra and the 
bottom panel shows the spectra scaled by Eq. ([3]). Al- 
though the error of the scaling relation is relatively large 
at high frequencies (> 10~^^ Hz), it is still reasonable 
at low frequencies where most of the energy is emitted. 
Thus, we use Eq. ([31) as a reasonable scaling relation for 
both equal- and unequal-mass mergers. 

The energy density of the GWB at redshift Zi can be 



III. RESULTS 

We obtain the spectrum of the GWB by using the 
method described in the previous section. First we 
show the dependence of the result with respect to the 
upper-cutofF mass Mcut of halos. Fig. |4] represents 
the GWB spectra for several mass cutoffs, IO^^/i^^Mq, 
IO^^H-^Mq, 1O"/i-1M0 and IO^^H-^Mq. The emitted 
energy reaches Qqw ^ 5 x 10"^" for Afcut < IO^^H-^Mq 
and the spectral shape is very similar to that of a single 
merger shown in Fig. [2] The difference in the spectra 
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FIG. 5: Contributions of several redshift ranges to the GWB 
spectrum. 
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FIG. 6: The total number of merger at each redshift. The 
assumed volume is (6 /i^^Gyr)^ correspoding to hubble scale. 



FIG. 7: The number of mergers at each redshift with respect 
to mass that is higher mass of two merger haloes. The as- 
sumed volume is (6 /i~^Gyr)'' correspoding to hubble scale. 
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FIG. 8: The spectrum by using ST(solid) and EPS(dashed). 



is very small for the cases with Mcut = l0^^^h~^AlQ and 
10^^h~^MQ. This can be easily understood by the shape 
of the mass function in Fig. [1] where the number of halos 
steeply decreases above 10^^^h~-^MQ. On the other hand, 
the difference is about two orders of magnitude between 
the cases with Mcut = IO^Mq and IO^^/i-^Mq, which 
implies that most of GWB energy is contributed from 
mergers of massive halos. Because the GWB spectrum 
is saturated at Afcut = IO^^/i-^Mq, we take lO^^h'^MQ 
as a cutoff mass hereafter. 

Next we investigate the relative contributions of red- 
shift ranges to the GWB spectrum. In Fig. [5l the contri- 
butions from several redshift bins are shown. The widths 
of the bins are determined so that each bin corresponds 
to the same cosmic time interval of 2 Gyr/h. Each spec- 
trum in Fig. [5] represents the GWB emitted by mergers 
in each redshift range. Due to the redshift of gravita- 
tional wave, contributions from larger z bins have peaks 
at lower frequencies. We find 46% of the total energy 
of the gravitational waves comes from the redshift inter- 



val < z < 0.19, while the gravitational waves from 
the 0.19 < z < 0.41 and 0.41 < z < 0.79 contribute to 
the GWB spectrum by about 26% and 15%, respectively. 
Thus, about 90% of the GWB comes from the redshift 
interval < z < 0.79. 

To understand this behavior, we show the merger rate 
below. Fig. [6] represents the total number of halo merg- 
ers for each redshift bin of Fig. [S] The total number has 
a peak at 0.79 < z < 1.39 and decreases toward lower 
z. However, as can be seen in Fig. [7] which represents 
the merger rate as a function of mass, mergers of mas- 
sive halos takes place more frequent for lower z. This 
would be natural in the hierarchical scenario of struc- 
ture formation. Because the energy of emitted gravia- 
tional waves increases substantially with the halo masses 
(see Eq. ([3])), the dominant contribution to the resul- 
tant spectrum of GWB comes from the lowest-redshift 
bin although the total number of mergers is relatively 
small. Here we remark that we have checked that the 
contribution from z > 2.2 is negligible as expected. 

Finally, we compare the spectra of the GWB calculated 



5 



using EPS and ST mass functions in Fig. |4l Here we took 
the cutoff mass Afcut = The total energy of 

the GWB for ST case is larger than that for EPS case 
by a factor of two and reaches il,GW ^ 10~^^. This is 
because massive halos are more abundant for ST than 
EPS. 



IV. SUMMARY AND DISCUSSION 

In this paper, we calculated the spectrum of grav- 
itational wave background originating from dark halo 
mergers by a quasi- analytic method. First, we con- 
structed merger histories by Monte-Carlo realizations for 
a mass function (Extended Prcss-Schechter formalism 
or the Shcth & Tormen formalism). Then we summed 
up the energy spectra from halo mergers following the 
merger history, using the result of jl2| which provides 
the gravitational- wave spectrum from a single merger as 
a function of the halo masses. We found that the energy 
density reaches Haw ~ 5 x 10~^° and that the domi- 
nant contribution comes from mergers of massive halos 
10^4 - IO^^H-^Mq at relatively low redshifts z < 0.19. 
We gave an interpretation of the relative importance of 
ranges of halo masses and redshifts showing the merger 
rates as a function of mass and rcdshift. We also com- 
pared the GWB spectra obtained by EPS to ST mass 
functions and found that the latter case has larger en- 
ergy by a factor of two. 

Finally, we discuss the observability of the GWB from 



halo mergers. Stochastic gravitational wave background 
converts the E-mode polarization of CMB into the B- 
mode through gravitational lensing between the observer 
and the last scattering. Thus observation of B-mode 
would be useful to probe GWs after last scattering. The 
energy of the GWB is at most ilcw ^ 10^^^ according to 
our calculations. This corresponds to inflationary grav- 
itational waves with the tensor-to-scalar ratio r ~ 10~^ 
0. It may be possible to detect them through B-mode 
polarization of CMB if the tensor-to-scalar ratio is lower 
than above the value. In this case, our result would be 
useful to probe the process of structure formation. 
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